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Comment on "Solidification of a Supercooled Liq- 
uid in a Narrow Channel" 

In a recent paper pi, Sabouri-Ghomi, et al. report 
on simulations of solidification in a channel geometry us- 
ing a phase- field method. They conducted two series of 
simulations, one with very small surface tension (relative 
to the channel width), and the other at larger surface 
tension. This system has been studied extensively both 
analytically and numerically j|, ||, §, |[ ||, 0| The au- 
thors do not relate to the more recent of these studies 
P, 0, 0, M which in fact explain many of the results the 
authors claim as surprising and can correct some misin- 
terpretations presented. 

Let us address the large surface tension case first. The 
authors find a transition from a tip widening instabil- 
ity for low undercooling, A ^ 0.5, to a stable finger for 
A > 0.5. The tip widening instability, as noted by Q] 
is a result of the lack of steady-state solutions for suf- 
ficiently low A. This threshold in A is a function of 
surface tension dp and anisotropy e. Thus, for example, 
at do = 0.01, (in units where the channel width L = 2) 
there are no steady-state solutions below A w 0.62 for 
zero anisotropy and below A w 0.60 at e = 0.1 [g. (Note 
that we define e as the anisotropy of the surface tension 
and not of the surface energy, as in Ref. [n|, which is a 
factor of 15 smaller for 4- fold anisotropy.) This thresh- 
old decreases with do , going to as do ^ and also with 
anisotropy. Above this threshold, there exist a bilater- 
ally symmetric finger, which is the channel analogue of 
the free dendrite. We have extended the calculation Q of 
the threshold to the parameters of the large surface ten- 
sion case of Ref. ||l[, namely (in our units) do = 0.00563 
and e = 0.75. The results are presented in Fig. 1, where 
we show the Peclet number, (dimensionless velocity), p, 
as a function of A. We see that the transition from the 
unstable Saffman- Taylor branch to the dendritic branch 
occurs at A = 0.49, in good agreement with the onset of 
stable dendritic growth reported in Ref. ||l|. The stable 
pattern seen above the threshold is thus nothing but the 
well-known channel analogue of the dendrite M, p|, and 
not a new type of solution. For example, we show in Fig. 
2 the dendrite for A — 0.55, which compares well with 
the pattern observed by the authors at this undercool- 
ing (see their Fig. 4). Also, this discussion should make 
clear that, contrary to the authors' speculations, there is 
nothing special happening at A = 0.5, which is of course 
the distinguished value in the Saffman- Taylor problem. 
We should also point out that at large A, there also ex- 
ist nonsymmetric fingers P, |l 0, 01; which have been 
termed parity-breaking dendrites, or doublons. The au- 
thors appear to have imposed reflection symmetry about 
the midline of their channel, which of course precludes 
seeing such solutions. 

For small do, the authors only find fingers which slow 
down in time, and widen correspondingly. The suprise 
here is that the authors did not see stable translating 
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FIG. 1: Peclet number vs. supersaturation A for symmetric 
solutions in a channel at do ~ 0.00563 and e — 0.75 
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FIG. 2: Dendrite for A — 0.55 for parameters as in Fig. 1 


dendrites. However, these simulations were performed at 
very small A, so it is possible that again the authors are 
below the threshold for dendrites. We in fact suspect that 
the effective do for this series of runs is in fact much larger 
than the nominal value do = 1.6 x 10~^ quoted, since it 
is extremely difficult to adequately resolve the interface 
over such small scales, even with adaptive gridding. 
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